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I. INTRODUCTION 

Numerical methods for ordinary differential equations are First the interval [a, b] is divided into an equal parts, x t = a 

methods used to find numerical approximations to the +ih, (i=0,1, 2 ...n], the step is h=x i+1 -Xi. 

solutions of ordinary differential equation .Their use is 
also known as "numerical integration”, although this term 
is sometimes taken to mean the computation of integrals. 


Many differential equation cannot be solved using 
symbolic computation. For practical purposes, however- 
such as in engineering-a numeric approximation to the 
solution is often sufficient. The algorithms studied here 
can be used to compute such an approximation. An 
alternative method is to use techniques from calculus to 
obtain a series expansion of the solution. 

Ordinary differential equations occur in many scientific 
disciplines, for instance in physics, chemistry, biology, and 
economics. In addition, some methods in numerical partial 
differential equation convert the partial differential 
equation into an ordinary differential equation, which 
must then be solved. 


A general form of the type of problem we consider is 
^=f(x,y ],x£[a,b] 

Y (a] = y 0 


Then to solve the function y(x] in a series of discrete 

equidistant node x 0 <x 1 <x 3 < . <x n to get approximate 

values y 0 <y 1 <y 2 <y3 < . <y n - 

II. PREDICTOR-CORRECTOR METHODS: 

The methods which we have discussed so far are called 
single-step methods because they use only the information 
from the last step computed. The methods of Milne’s 
predictor-corrector, Adams-Bash forth Predictor corrector 
formulae are multi-step methods. 

In solving the equation — =f(x, y], y (x 0 ] =y 0 we used 


Euler’s formula 

yi+i=yi+h/'(Xi,yi],i=o, l, 2 .(i] 

We improved this value by Improved Euler method 
yi+i=yi+l h [f(x;,y;] +f [X i+ 1 ,y i+1 )] .(2] 


In the equation (2], to get the value of y i+1 we require y i+1 
on the R.FI.S. To overcome this difficulty, we us we predict 
a value of y i+1 from the rough formula (1] and use in (2] to 
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correct the value. Every time, we improve using (2], Hence EULER’S METHOD 

equation (1] Euler’s formula is a predictor and [2] is a Euler’s method is the simplest one-step method and has a 

corrector. A predictor formula is used to predict the value limited application because of its low accuracy. However, 

of y at x i+1 and a corrector formula is used to corrector the it is discussed here as it serves as a starting point for all 
error and to improve that value ofy i+1 . other advanced methods. 

2.1. Milne’s Predictor Corrector Formula Given the differential equation 

The Milne-Simpson method is a Predictor method. It uses —=f(x, y) 

a Milne formula as a Predictor and the popular Simpson’s dx 

formula as a Corrector. These formulae are based on the . ... , , , 

. , , , . , , Given the initial condition y fx 0 ) = y 0 

fundamental theorem of calculus. 

Y (Xj+iJ = y 0,0 + J Xi+1 f(x,y)dx .[3] dy £ , 

r i+u j c jj j Xj cj Wehave^=/(x 0 ,y 0 ) 

When j=i-3, the equation .becomes an open integration And therefore Y[x ) =y (je )+ [x . Xq] /(*„, y Q ) 
formula and produces the Milne’s formulae 

yt+i, P =yt-3+ Y( 2 /t- 2 ‘/t-i +2 /t] [4] Then the value of y(x] at x=x 1 is given by 

Y fo) = y [x^y (x 1 -x 0 ) f(x 0 , y 0 

Similarly, When j=i-l, equation becomes a closed form 

integration and produces the two-segment Simpson’s Letting h=x 1 -x 0 , we obtain y 1 =y 0 +h/(x 0 , y 0 ) 
formula 

Vt+i ,c=yi+i (/i-i+T/i+Zi+iJ (5) Similarly y(x) at x=x 2 is given by 

3 y 2 =yi + h/(^i-yi) 

Milne’s formula is used to‘predict’the value of y, +1 which . , , . . 

. , , , , . r In general, we obtain a recursive relation as 

is then used to calculate G +1 , r , , „ , _ 

f _ frv „ y n+1 =y„+h/(x n ,y n );n=0,1,2. 

Ji+l-J f x i+l< Xi+lJ 

_ Improved Euler method 

Then, equation (6) is used to correct the predicted value . . ., . . ., ,. . ., , . , ., 

r „ . , , Let the tangent at(x n , y n ) to the curve be p n A. In the 

ofy, +1 . The Process is then repeated . . ,, 6 , , v . „ , , j 

r \ , r . „ , . , r , . , , . interval!*,,,x,), by previous Eulers method, we 

for the next value of i. Each stage involves four basic calculations, namely, , ' , ^ ^ . 

„ „ . ° approxfrnate the curve by the tangent p„A. 

1. Prediction of y i+1 y 6 

3. correction ofy- + + i Ti® =%+¥(*„, y 0 ) Where y^=M lQl 

4. improved value of /y +1 [for use in next stage ] 

y Q 1 [x 1 ,y 1 (1) ). Let Q 1 C be line at Q t whose slope is f 

It is also possible to use the corrector formula repeatedly t Xl> ^ Now take the average of the slopes at P 0 and Q 1 
to refine the estimate of y i+1 before moving on the next Le ' 
stage. -[/(*o- y 0 ]+f C x lt y/ 15 )] 

2.2. Adams-Bash forth-Moulton Method Now draw line p 0 D through P 0 ( x 0 ,y 0 ] with this as the 

Another popular fourth-order predictor-corrector method slope. 

is the Adams-Bash forth-Moulton Method multistep Thatis,y-y 0 = ~[f(x 0l y 0 ) +f ( x^y^)] (x-x 0 ) 

method. The predictor formula is known as Adams- Bash 2 

forth predictor and is given by _ . 

h os This j ine in tersects x=x, at 

ym = yi+^(55/ r 59/ i _ 1+ 37/ i _ 2 -9/ i _ 3 ) .(6) y 1= y 0+ i/r[/(x 0 ,y 0 ) +f(x 1 ,y 1 «J] 

1 

The corrector formula is known as Adams-Moulton yi-yo + “/i|/(^o<yo) + /(*i<yo + h/(Ao<yo))] 

Corrector and is given by 

yi+i = y t +£ ifi-z-Sfi* .(7) Writing generally, 

yn+i=yn + -h[f(x n ,y n ]+f(x n + h,y n +h ( , ]]] 

This pair of equations can be implemented using the 

procedure described for Milne-Simpson method. RUNGE-KUTTA METHODS 

Runge-Kutta method reduces data requirements to reach 

TAYLOR SERIES METHOD the same precision without higher derivative calculation. 

In mathematics, a Taylor Series is a representation of a 

function as an infinite sum of terms that are calculated The numerical solutions of the numerical equation 

from the values of the functions derivatives at a single —=f( X; y) 

point 

The numerical solution of the equation Given the initial condition y (x„) = y 0 the fourth order 

—=f(x y) Runge-Kutta method algorithm is mostly used in problems 

dx unless otherwise mentioned. 

Given the initial condition y (x 0 ) = To k^h f (x, y) 

yi=y (*i)=yo+^yo'+^-yo" + . k 2 =hf(x + -h,y + -k 1 ) 
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k 3 =hf(x + \h,y + jfe 2 ) 
k 4 =hf(x +\h,y + \k 3 ) 
Ay=^( fe i + 2fc 2 + 2fc3 + M 


EXAMPLE 

1. Determine the value of y (0.4) using Milne’s method 
given y'=x+y, y (0) = 1 use Taylor series to get the values 
ofy (0.1), y (0.2) andy (0.3). 


Solution. Here x o =0, y 0 =1, Xj =0.1, x 2 =0.2, x 3 = 0.3, x 4 =0.4 


y — x + y 
y'= 1 +y 

y'"=y" 


y 0 =x 0 +y 0 =0+1 =1 
yo"=i + yo=i + i =2 

To = To = 2 


Y (0.1) = 1+ (0.1) (1) +0.01+0.0003 =1.1103 

Y (0.2) = 1.1103+0.1210+0.0111+0.0003 =1.2427 

Y (0.3) =1.2427+0.1443+0.0122+0.0004 =1.3996 
Now, knowing y 0 ,yi, y 2 ,y 3 we will findy 4 . 


By Milne’s predictor formula 

y 4 , p =l+0.1333[2.4206-1.4427+3.3992] 
yi'=Xi+yi =0.1+1.1103 =1.2103 
y 2 '=x2+y 2 =0.2+1.2427 =1.4427 
y 3 ,=:e 3+y3 -0.3+1.3996 =1.6996 
y 4p =1.5835 

y 4 , p ’=x 4 +y 4 =0.4+1.5835 =1.9835 

Milne’s Corrector formula 

y 4c =l.2427+0. 0333(1.4427+6.7984+1.9835] =y (0.4) 
=1.5832 

Adam’s predictor formula 

y 4p = 1.3996+0.0042(93.478-85.1193+44.781-9] =1.5849 
y 4 , P '=x 4 +y 4 =0.4+1.5849 =1.9849 

Adam’s corrector formula 

y 4c = 1.3996+0.0042(17.8641+32.2924-7.2135+1.2103] 
=y (0.4) =1.5851 2. Determine the value of y (0.4) using 
Milne’s method given y'=x+y, y (0) =1 use Euler’s 
method to get the values of y (0.1), y (0.2) and y (0.3). 


yn+l=yn + h/On.yn] 

y 4 = 1+ (0.1) /(0,1)= Y (0.1) =1.1 

y 2 = 1.1+ (0.2) /((0.1), (1.1)= Y (0.2) =1.22 

y 3 = 1.22+ (0.1) /(0.2 4- 1.22)= Y (0.3) =1.362 


By Milne’s predictor formula 

y 4p =1+0.1333 [2.4-1.42+3.324] =1.5737 

yi'=*i+yi =0.i+i.i =1.2 

y 2 '=x 2 +y 2 =0.2+1.22 =1.42 
y 3 '=x 3 +y 3 =0.3+1.362=1.662 
y' =1.9737 

s 4 ,p 


Milne’s corrector formula 

y 4 c =l.22+0.0333(1.42+6.648+1.9737] = y (0.4) =1.5544 


By Adam’s predictor formula 

y 4 p =1.362+0.0042(91.41-83.78+44.4-9] =1.5427 
y'. =1.9427 

J 4 ,p 


Adam’s corrector formula 

y 4c =1.362+0.0042 [17.4843+31.578-7.1+1.2] =y (0.4) 
= 1.5433 

3. Determine the value of y (0.4) using Milne’s method 
given y'=x+y, y (0) =1 use Improved Euler’s method to 
get the values of y (0.1), y (0.2) and y (0.3). 

yn+i=yn + lk[f(x n ,y n )+f(x n + h,y„+h/(x n ,y n )] 

yr=l+ ko.l) [/(o, 1) + /(0+0.1, 1+0.1/(0, 1)] = Y (0.1) 
= 1.11 

y 2 = l.ll+ ^(0.1) [/(0.1, 1.11) + /(0.1+0.1, 1.11+0.1/(0.1, 
1.11)] = Y (0.2) =1.2421 

y 3 = 1.2421+ ^(0.1) [/(0.2, 1.2421) + /(0.2+0.1, 

1.2421+0.1 /(0.2,1.2421)] = Y (0.3) =1.3985 

By Milne’s predictor formula 

y 4 p=l+0.1333[2.42-1.4421+3.397] 
y 1 '=x 1 +y 1 =0.1+1.11 =1.21 
y 2 '=x2+y 2 =0.2+1.2421 =1.4421 
y 3 ’=x 3 +y 3 =0.3+1.3985 =1.6985 
y 4 p =1.5832 

y 4 ,p'=x 4 +y 4 =0.4+1.5835 =1.9832 

Milne’s Corrector formula 

y 4c =l.2421+0. 0333(1.4421+6.794+1.9832] =y (0.4) 

= 1.5827 

Adam’s predictor formula 

y 4p = 1.3985+0.0042(93.4175-85.0839+44.77-9] =1.5837 
y 4 , P '=x 4 +y 4 =0.4+1.5849 =1.9837 

Adam’s corrector formula 

y 4c = 1.3985+0.0042(17.8533+32.2715-7.2105+1.21] =y 
(0.4) =1.5838 

4. Determine the value of y (0.4) using Milne’s method 
given y'=x+y, y (0) =1 use Runge-Kutta method to get 
the values of y (0.1), y (0.2) and y (0.3). 

k t = (0.1) (0+1) =0.1 
k 2 = (0.1) (0.05+1.05) =0.11 
k 3 = (0.1) (0.05+1.055) =0.1105 
k 4 = (0.1) (0.1+1.1105) =0.12105 
Ay=0.11034 

Y (0.1) =1.11034 

k 4 = (0.1) (0.1+1.1103) =0.1210 
k 2 = (0.1) (0.05+1.05) =0.1321 
k 3 = (0.1) (0.05+1.055) =0.13264 
k 4 = (0.1) (0.2+1.1105) =0.1443 

Y (0.2) =1.2428 

k 4 = (0.1) (0.2+1.2428) =0.1443 
k 2 — (0.1) (0.25+1.3149) =0.1565 
k 3 = (0.1) (0.25+1.3211) =0.1571 
k 4 — (0.1) (0.4, +1.3999) =0.1799 

Y (0.3) =1.4014 
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By Milne’s predictor formula Milne’s Corrector formula 

y 4p =l+0.1333 [2.4207-1.4428+3.4028] y 4c =1.2428+0.0333[1.4428+6.8056+1.9839] = y [0.4] 

yi'=*i+yi =0.1+1.11034 =1.2104 =1.5835 

yz = x 2+y 2 =0.2+1.2428 =1.4428 Adam’s predictor formula 

y 3 '=x 3 +y 3 =0.3+1.4014 =1.7014 y 4 _ p = 1.4014+0.0042[93.577-85.1252+44.78258-9] 


y 4p =1.5839 

y 4iP '=x 4 +y 4 =0.4+1.5835 =1.9839 


= 1.5872 

y 4 , p '=x 4 +y 4 =0.4+1.5872 =1.9872 


Adam’s corrector formula 


y 4 c = 1.4014+0.0042[17.8848+32.3266-7.214+1.21034] =y [0.4] =1.5871. 



X 

Taylor’s series 

Euler’s method 

Improved Euler's method 

Runge-Kutta method 

Exact Solution 

0 

1 

1 

1 

1 

1 

0.1 

1.11034 

1.1 

1.11 

1.11034 

1.1103 

0.2 

1.2427 

1.22 

1.2421 

1.2428 

1.2428 

0.3 

1.3996 

1.362 

1.3985 

1.4014 

1.3997 

n a 

Milne’s 

1.5849 

1.5544 

1.5827 

1.5835 

1.5836 

U.4 

Adams 

1.5851 

1.5433 

1.5838 

1.5871 

1.5836 


Error Value 



X 

Taylor’s series 

Euler’s method 

Improved Euler’s method 

Runge-Kutta method 

0 

0 

0 

0 

0 

Ss 

: 0.1 

-0.00004 

0.0103 

0.0003 

-0.00004 

0.2 

0.0001 

0.0228 

0.0007 

0 

0.3 

0.0001 

0.0377 

0.0012 

-0.0017 

n a 

Milne’s 

-0.0013 

0.0292 

0.0009 

0.0001 

U.4 

Adams 

-0.0015 

0.0403 

-0.0002 

-0.0035 


CONCLUSION 

In this paper, we compared the various numerical method 
with the Adam’s predictor and corrector and Milne’s 
predictor and corrector. The problem of Runge-Kutta 
method has the minimum error value. So the Runge- 
Kutta method problem is best to approximate and exact 
solution. 
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